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1 Introduction 

Following leading examples of mathematical physics, abstract linear evolutionary problems of 
the form 

doV + AU = f on R (1) 

come into focus. Here do denotes the operator of time-differentiation, established as a normal, 
boundedly invertible operator in a suitable Hilbert space (see |13|). and we assume that A is 
a closed, densely defined operator, such that ^4 + Aq and A* + Aq are both maximal accretive 
for some Aq £ M in some Hilbert space setting. To determine V and U from the given data / 
the equation ([1]) needs of course to be completed, here by an additional rule of the form 

V = MU (2) 

frequently referred to as a "material law". Here Ai is given in terms of an operator-valued 
function calculus associated with the (inverse) time-derivative via a bounded-operator-valued 
and analytic function M. The theory is set up in |10| for the case A skew-selfadjoint and in 
|12| [T6] for the general case resulting in well-posedness of the evolutionary problem (we omit 
the closure bar over the operator sum) 

{doM {d^^) +A)U = fonR. (3) 

In order for this to describe a process evolving in time we not just want the existence of a 
solution operator [OqM (^d^ ) -|- j4) but we would also expect a causality property to hold. 
Causality, as the property that solutions vanish as long as data are zero, can be conveniently 
encoded as 

X]_,.[ (mo) {doM {d^') +A)-' = xj_,„, (mo) {doM {d^^) + Ay' X]_,„[ (^o) (4) 

for all a € M. Here X, ^ , (m-o) denotes the temporal cut-off operator 

(Xj_,.[ (mo) ^) (t) = X,_,„[ (t) ^ (t) for t G M. 

The symbol niQ serves as a reminder for 'multiplication-by-argument' with respect to the time 
parameter. Since the systems considered are time-shift invariant it would suffice to require 
@ just for one a € M - say a = 0. 

In this approach well-posedness is warranted by the eventual (real) strict positive definiteness 
of doM (^Oq ) for all sufficiently large g G ]0, cxd[. In this paper we want to consider material 
laws of the specific form 

M (Sq-I) = Mo + ^ 5o-"M, + d^Hli + d^^Ah {do') , (5) 

where H C ]0, 1[ finite. Ma continuous and selfadjoint for a S H. The terms Oq" are the 
fractional integrals of the title, which can be defined in terms of the said function calculus. Such 
material laws have been of interest in many applications for a long time and although the strict 
positive definiteness requirement is exhaustive, it is not clear what manageable conditions on 
the coefficient operators Mq, M^, a € H, would result in a corresponding positivity estimate, 



2 Functional analytic framework 



see ()13p below. It is the purpose of this paper to give such conditions. For this we first 
need to estabhsh the said function calculus associated with the time derivative do allowing 
in particular for a proper definition of the fractional integral Oq'^, a G ]0, 1[. This is done in 
the Section 12.11 Although there exists a vast literature on the topic of fractional derivatives 
(see for example the monographs |14l H] and the references therein), the article is largely 
self-contained, so that there is no much need to refer to classical statements on fractional 
derivatives to understand the results in this work. To the best of the authors' knowledge it 
has not been widely noted that do can be established as a normal operator, see [T^, which then 
comes with its own standard function calculus via the spectral theorem for such operators. 
Indeed, it turns out that a unitary variant of the Fourier-Laplace transformation yields a 
spectral representation for 3xndo = 2i(^0 ~ "^o) ^^^ 9^c9o is just multiplication by a number. 
This simplifies matters considerably and since we are staying in a Hilbert space setting there 
is indeed no need to utilize more intricate results of fractional calculus in other spaces. This 
approach was already successfully used to study homogenization in fractional elasticity, |17| . 
After establishing the function calculus we shall characterize in Section [3] a class of material 
laws for which ()13p can be shown to hold. This is the main result of the paper. Compared 
to the results in |17l Theorems 2.1 and 2.2] and |9] we will find that our more in-depth study 
here implies less restrictive assumptions on the class of material laws. The final Section H] 
is devoted to illustrating the general results with applications to a fractional Fokker-Planck 
equation (Section 14. ip and a fractional Kelvin- Voigt type model from solid mechanics (Section 
14. 2p . Section [J] is concluded with some remarks of how to deal with initial value problems 
f Section [Op . 



2 Functional analytic framework 

2.1 Fractional Calculus and Operator- Valued Functions of Time 
Differentiation 

As indicated in the introduction, we start by establishing time differentiation Oq as a nor- 
mal operator. We consider the weighted f/^- valued L^-type space Hgfi{^,H), generated 
by completion of Coq{^,H), the space of smooth i7- valued functions with compact sup- 
port, with respect to the inner product ( • | ■ )poo given by the weighted integral (y'l^) ^ 
Jjg {(f (t) {ip {t))ff exp {—2Qt) dt. The associated norm will be denoted by | • | „ g O' '^^^ multi- 
phcation operator Coo (IR, H) C Hq^ (M, H) -^ Coo (I^, H) C Hq^q (M, H) = L^ (M, H) , ip ^ 
(t >-^ exp { — gt) if (t)) , clearly has a unitary continuous extension, which we shall denote by 
exp ( — qitiq). Its inverse (adjoint) is given by exp {gmo). By taking the closure of the operator 

Coo{R,H) C H,,o{^,H) ^ H,,o{^,H) 
(p^cf)', 

the time-derivative do can be established as a normal operator on Hq^q (M, H) with real part g 
and so with ^ {do — g) as imaginary part. The domain of do can be characterized by functions 
belonging to Hgfi{M.,H), whose weak derivatives also lie in Hg^oi^yH). For f? € M \ {0} we 
have the bounded invertibility of do- The inverse of the normal operator do is bounded and 



2.1 Fractional Calculus and Operator-Valued Functions of Time Differentiation 
can be described by 



[- Jt fields if^<0 



for every / £ -ff^.o (Mi H) and almost every t € M as a Bochner integral. Henceforth we shall 
focus on the case q S ]0,oo[, which is associated with (forward) causality, see e.g. equation 
dU. The Fourier-Laplace transform C^ := J-"exp( — ^mo) : HgfiCRjH) — )• L'^{W,H), given 
as a composition of the (temporal) Fourier transform J-" and the unitary weight operator 
exp { — gmo), is a spectral representation associated with Oq. It is 

do = £* (imo + £') Cg, 

where niQ denotes the multiplication-by-argument operator given as the closure of 

Coo (K, H) C L^ (M, H) -^ L^ (M, H) 
ip I-)- mQip 

with 

{rriQip) (X) ■.= X{p{X) in H 

for every A € M. This observation allows us to consistently define an operator function 
calculus associated with do in a standard way, [2]. Clearly, we can even extend this calculus 
to operator-valued functions by letting 



*'(*') ■="''" (t^) "'■ 



Here the linear operator M ( ^^ , j : L^ (M, H) — t- L^ (M, H) is determined uniquely via 

for every A € M, 93 G Coo (IR) H) by an operator-valued function M. For a material law 
the operator-valued function M needs to be bounded and an analytic function z 1— )• M (z) 
in an open ball -Be {i", f) with some positive radius r centered at r. This is not an artificial 
assumption, rather a necessary constraint enforced by the requirement of causality, see [10] or 
|15| Theorem 9.1]. 

In terms of the associated operator- valued function calculus we also know what 

Sq-", ae[0,l[, 
meano With this we define for 7 G M 

5?:=4^^9j-r^^ (6) 



^It should be noted that Oq " is largely independent of the particular choice oi g € ]0, cxd[. Indeed, since 



2 Functional analytic framework 

as a natural generalization of differentiation to arbitrary real orders. Here [a] denotes the 
smallest integer greater or equal to a. The fact that (Oq) ^ is a family of commuting operators 
appears to be useful in applications, see |6] and the quoted literature. In contrast, the varianin 

aDl ■■= 5JX,_[ ("^O) = 4"^ (9j"^"^X,_i (-0) 

with an appropriate choice of domain, known as Riemann-Liouville fractional derivative, a E M 
a parameter, |14| H], lacks this property. This is also true for the frequently used alternative 
fractional derivative, the Caputo fractional derivative, [TJIU]: 

with a suitable domain. With these fractional derivatives being mere variants of limited 
usefulness in our context, we shall utilize only the above spectral definition ([6]) for our fractional 
differentiation/integration. 

Let us inspect more closely some properties of 5q for a S ]0, 1[. Denoting by || • Hj,^ the 
operator norm for operators in the space X we record our first lemma. 

Lemma 2.1. For a G [0, 1] we have 

^td^>Q'' (7) 

and for a G [0, oo[ 

ll^cr"Lg,o(R,H) - s~"- (8) 

Proof. To see ([7]) it suffices - by the function calculus of do - to consider IHe (it + g)" for t G M. 
We have 

me ((it + gf) = lit + ^|" cos (a/3) 

a/2 



Q 



1 + (tan /3) j cos (a/3) 



with 



(cos p) 



/3:=arg(it+^)G]-7r/2,^/2[. 



we have for ^p € Coo (K, H) 



1 



and 






From this convolution integral representation we can also read off that d^ " is causal. 
■^In the limit a — !> —00 the spectral fractional derivative is formally recovered: 

^7 _ r)7 



2.1 Fractional Calculus and Operator-Valued Functions of Time Differentiation 

The values of ^"^^^fi are either unboundeqfl as /3 — > ib^ (a g]0, 1[), or 1 (if a € {0, 1}). In 
the first case we shall try to find a minimum. The derivative 

cos (a/3) \ I ^ a (cos/3)° sin (a/3) — cos (a/3) a (cos (/3))"~ sin (/3) \ 

(cos/3) J y (cos/3)^" J 

cos /3 sin (a/3) — cos (q/3) sin (/3) 
(cos/3)"+^ 

sin((a-l)/3)' 



/3 I-)- -a 

/3 IH> -Q . -, 

(cos/3) 

vanishes only for /3 = 0. Due to the behavior at the limits this must be the point where a 
minimum occurs. For /3 = we have, however, 

cos (a/3) 
(cos/3)" 

which proves ([7]). Estimate ([8]) now follows for a G [0, 1] since 

l^l.,o,o 1^0 t^l,,o,o > ^^ (f^l^o f^),,o,o > £-" (f^|f^).,o,o 

and so 

|5o"f/|,,o,o>^"|f^l.AO 
for allU € D (d^), which implies ([S]) for a G [0, 1]. Since 

with [aj denoting the largest integer less or equal to a, ([8]) follows for arbitrary a G [0, oo[ by 
the submultiplicativity of the operator norm and the fact that Wd^ \\ < g~^. D 

For our purposes we mostly need to be concerned with a G ]0, 1[ leading to 

d^ = dod^-\ae]o,i[. 

Lemma 2.2. Let I be a compact interval in ]0, 1[. Then the function 

with unbounded selfadjoint operators IHe (9q ) in Hg^Q{M.,H) as values, is monotonically in- 
creasing in the sense that there exists Qq > such that for all Q > Qo o,nd for all a < (3 in I 
the estimate 

holds for all U & D{d^). 

Proof. We define 

^0 '■= exp ( — tan ( max/ — ) ) > 0. 



'This fact makes the cases a G {0, 1} and a G ]0, 1[ incomparable. 



2 Functional analytic framework 

Let t G M be fixed and g > Qo- We recall from ([9]) that for a G [0, 1] 

with 

/3:=arg(it + ^)G]-7r/2,7r/2[. 

Consider now 

" i-^' q""-, TWa I = I a 1-^ I t; ] cos (a/3) ( In - /3tan (a/3) 

(cos/3) / \ \cosp/ \ cosp 



Noting that 



m — p tan (ap) > m o tan a— 

cosp 2 \ 2 



> In rt tan max / — > 

- '^ 2 V 2/ 

for all a G / (note < max/ < 1) we see that for a < /3 in / 

IRe ((it + £))")< ^e ((it + ^)^) (10) 

for all t G M. Via the function calculus of do this yields the (dense) inclusion 

For U e D (d^) we also see from dlOD that 

{U\ D\t (d^) [7)^,0,0 = {^,U\ 9\t ((imo + qT) C,U),,, 

< {C,U\ 9.e ((i.^0 + ,r) C,U)^^^^^ = {U\ 9.e {d^) u) ^^^^^ . D 

2.2 Sobolev-chains 

In this subsection we recall the notion of Sobolev chains (see |10| Section 2.1]) associated with 
a closed, boundedly invertible operator on a Hilbert space. Throughout let H he a Hilbert 
space and C : D{C) Q H —^ H a densely defined closed linear operator with G q{C). 

Definition 2.3. For A: G Z we define the inner product {■\-)Hk(C) on D{C ) by 

{Ay)H,iC) ■■= {C''x\C^y)H (x,y G Hk{C)). 

Moreover, we define Hk{C) as the completion of D{C ) with respect to the inner product 
{■\-)Hk{C)- The sequence (-fffc(C'))fcez is called the Sobolev-chain associated with C. 

Proposition 2.4 ([TO', Theorem 2.1.6] ). For each j,k £ Z with j < k we have Hk{C) ^^> 
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Hj(C), where the embedding is continuous and dense. Moreover, the operator 

X I—)- Cx 

has a unitary extension, which will he again denoted by C. 

Example 2.5. For £> > the operator do : D{do) C Hgfi{R,H) -^ Hg^o{M,H) is densely 
defined and closed. Moreover, since 9^e9o = Q^ we have G Q{do). Thus, we can define the 
Sobolev-chain associated with the time-derivative, which will be denoted by 

Then the Dirac-distribution 6 lies in Hg^^i{R,H) with Oq 6 = X[o,oo[ ^ Hg^CK,!!). 

Remark 2.6. We can extend the operator C : D{C) Q H ^ H to the Hilbert space Hg^Q(M., H) 
for each ^ > in the canonical way, i.e., we set {Cu) (t) = C {u{t)) for every u € ifg^o(I^) H) 
taking values in the domain of C and almost every t € M. Then, obviously the operators Oq 
and C commute. Thus, for each k,j G Z the operators can be realized as unitary operators 

5o : if,,fc+i(]R,i/,-(C)) ^ /7,,fc(M,if,_i(C)) 

and 

C : H,^k{%H,+i{C)) ^ H,^k{%H,{C)). 

Moreover, for all M E L{Hgfi{M.,H)) being such that d^ M = MOq there is a unique 
continuous extension to an operator on Hgi^(E.,H) for all A: G Z. 

3 A Class of Fractional Material Laws 

We are now able to rigorously establish the specific form of the material law operators we wish 
to consider, namely 

M (Sq-I) =Mo + J2 ^o""^« + do^Mi + d^^h (Sq-I) . (11) 

aen 

Here 11 C]0, 1[ is a finite set, Mq, for a € {0} UlIU {1} are bounded linear operators in H and 
M2 is a bounded operator valued function, analytic in a ball Bq (r, r) for some r G ]0, cx)[ with 
r > 2^ , compare \W\ [T2] . Then the well-posedness of the corresponding evolutionary problem 
([3]) can be shown with the help of the main result presented in [7]. 

Theorem 3.1 ([7, Solution Theory] ). For r > let M : Bc{r,r) — > L{H) be a bounded 
analytic mapping satisfying 

\J /\ d\ez-^M{z)>c. (12) 

c>0 zeBc[r,r) 

Moreover, let A : D{A) Q H ^ H be a skew-self adjoint operator. Then, for each Q > 2^. the 
operator doM{dQ ) -\- A is invertible and the inverse 
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3 A Class of Fractional Material Laws 

is bounded and causal in the sense of Q). Moreover, for each k & Z the solution operator 
(5oM(9(^ ) + Aj can be established as a bounded linear operator on Hg^i.(W,H). 

In order to apply this well-posedness result, we need to require the positive definiteness con- 
straint (|12p for our material law of the form (lllll . that is for some constant cq € ]0, cxd[ and all 
sufficiently large g € ]0, oo[ we require 

IRe {U\doM {d^') U)^^^^^ > CO (t/|f/),,o,o (13) 

for all U (z D [Oq) by the unitarity of the Fourier-Laplace transformation. 

If (|13p holds for the case M2 [Oq ) = with a constant cq € ]0, oo[ and for all sufficiently large 
Q S ]0, oo[, we may consider a non- vanishing remainder term M2 (d^ ) as a perturbation by 
assuming 

limsup||M2 (5o"^)||„ ^(K,H) < CO, (14) 

yielding an estimate of the form (|13p with cq replaced by a positive constant cq < cq — 
limsup^.,^ ||M2 (^o"^) L,,o(R,iT)- 

Our aim is to give conditions in more explicit terms which warrant condition ()13p . 



For convenience we shall use a monotonically increasing enumeration (ao, • • • , oat) of 11. For 
the purpose of our following considerations we first record the following rather elementary 
observation. 

Lemma 3.2. Let l be the canonical (isometric) embedding of a closed subspace V ^ H into 
H. Then a* is the orthogonal projector onto V . Let k be the canonical embedding ofV into 
H then we have 

U* + KK* = 1. (15) 

Proof. It is 

X I— )• x 

an isometry with the inner product {■ \ ■ )y taken as just the restriction of the inner product 
{■ \ ■ )jj oi the Hilbert space H. We find 

{x\y)H = {i'^\y)H = (^k*y)y 

for all X G y and y ^ H. We read off that for y £V we have 

{x\y)y = {x\y)j^ = {ix\y)H = {A'*y)v 
for all X G F and so 

i^*y = y- 

Moreover, for y G V^ we read off 



= {x\i*y)y 
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for all X G V, i.e. 

L*y = 0. 

Consequently, we also have 

ii*y = y for y G V, 
u*y = for y G V-^. 

These properties, however, characterize the orthogonal projector onto V. Similarly for k and 
the projection theorem yields ()15p . D 



Remark 3.3. Equality (jlSp may be written in an intuitive block operator matrix notation as 

(.«)(f.)=L 

It may also be worth noting that l* i : V ^ V and k* k : V — )■ V are just the identities on 
V and V-^ , respectively. It is common practice to identify H and V © V^ , which makes 

the identity. For our purposes it appears helpful to avoid this identification. The mapping 
(|16p is obviously unitary, which allows us for example to study an equation of the form 

NU = F 



for a bounded linear operator N in H via the unitarily equivalent block operator matrix 
equation 

::)~("^)(::)-4::::::::)(:".)-(::". 



Note that 

is the inverse of (11611. 



Let Pq, Fq, Qq orthogonal projectors in H with 

Po + Fo + Qo = l 
and 

the corresponding factorizations in the above sense. We shall impose the following general 
assumption. 

Condition 3.4. Let the continuous linear operator Mi and the selfadjoint, continuous oper- 
ators Mq, Mq^, k E {0, . . . ,N}, be such that Mq, M^^ are commuting with Pqj Qo (and so 
also Fq). Moreover, let Po-^Ofc-Po) QoMaf.Qo > 0, A; € {0, . . . , A^}, and Mq non-negative and 
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3 A Class of Fractional Material Laws 

such that 

L*p^MQLpg, Lg^yit MiiQg, L*p^MaQLFo are strictly positive definite. (17) 

Theorem 3.5. Let Condition\3.4\ hold. Then the well-posedness requirement lilcl\) for 



M(9o"i) = A/o + Y, do'^Ma + Oq'Mi 
aen 



is satisfied. 



Proof. By the Fourier-Laplace transformation, it suffices to estimate £'Afo + X]A;=o '^" (s) ^a^ + 
D\t Ml uniformly with respect to t G M and all sufficiently large g € ]0, cx)[. Here, we have set 

Due to the assumed commutativity of Pq and M^ and the non-negativity of P^MaPo = 
PoMa = M^Pq and QoM^Qo = QqM^ = M^Qo, a G 11, we get with some generic constants 
CO) ci, C2 € ]0, oof using standard estimates of the form 2ab < etP' + e~^l?' that 



N 



qMo + Zk=o </>«, {q) M„, + D\t Ml 

= gPoMoPo + gFoMoFo + gQoMoQo + ELo '^"fe is) PoMa.Po 

+Po d\t MiPo + Fo 9^e MiFq + Qo 9^e MiQo 

+2 IHe Fo 9^e MiPq + 29\eQo 9\t MiPq + 2D\eQo D\c MiFq (18) 

> £.PoMoPo + E^=o </>«. (e) ^oM„,Fo 

+Po 5^e MiPo + Fo 9^e MiFq + Qo 9^e MiQo 

+2 IRe Fo IRe MiPo + 2 IHe Qo 9^e MiFo + 2 «He Qo 9^e MiFo 

> £»coFo + Y.k=o (t>ak {q) FoMa^Fo - C2F0 + CiQo 

holds for all sufficiently large ^ G ]0, oo[. In consequence we may reduce our considerations to 

N 



^(/)„, (e)FoM«,Fo-C2Fo. 



k=0 

Indeed, using (fTO|) we obtain 

Ylk=o(t>a,{e)l'*FoMa,LFo > 

> (l)ao (q) C3 + T.k=l (t>ak (q) l'*FoMa^l^Fo (19) 

for some C3, C4 > 0. Thus, from (|19p we get for all sufficiently large g G ]0, cxd[ 

N 



J24>akl'*Fo^<^k''Fo > y0ao {q) , 



fc=0 



where we use that 



"^-(^^ ^0 (,^oc) 



^ao 



(0) 
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by ([9]). Using that according to Lemma l2. II we have (pao (q) ^ Q^ "" ^^'^ inserting this into 
(pTSj) yields 



N 



qMo + ^ ,/.„, (^) M„, + $He Ml > ^cqPo + ( y^^""" - ci) Fq + C2Q0 

fe=0 

> min I t'Co , y ^i^^"" - ci , C2 1 

and so the desired strict positive definiteness holds indeed for all sufficiently large g G ]0, oo[. 

D 



We may summarize our findings in the following well-posedness theorem, which according to 
the above is just a particular case of the main result of [7] (see Theorem 13. ip . 

Theorem 3.6. Let M (8^ ^ be a material law of the form ^ satisfying Condition \3.4\ and 
( jj^p . Then for every sufficiently large g € M>o and every f G Hgi^. (M,H) there is a unique 
solution U G Hg^k 0^,H), /c G Z, satisfying equation ^. Moreover, the solution operator 

{doM {do^) +A)-^: H,^k (M, H) ^ i/,,,. (R, H) 

is continuous and causal in the sense that, for every a G M, if g & Hg^ {^^ H) vanishes as an 
H-valued generalized function on ] — 00, a[ then so does [OqM (SJ" ) + A) g. 

Remark 3.7. 

1. Causality implies that if 5 G -^^,0 (I^; H) vanishes on ] — 00, a [ then the fractional deriva- 
tives occurring in the material law are actually Riemann-Liouville type derivatives aD'j^'' 



t I 



A;G{0,...,A^}. 



2. We have chosen to discuss the case of A skew-selfadjoint in H. More general operators 
A could also be admitted by using other suitable boundary conditions, compare e.g. 
|1HI16). The operator A can even be generalized to a monotone relation, see |16) . 



We shall give a first illustrative example in which case the above theorem may be applied. 

Example 3.8. Let Hq,Hi,H2 be Hilbert spaces. For a,/3 g]0, 1[ with a < /? we consider the 
following operator matrices with respect to the Hilbert space Hq © Hi © H2: 



Mn := I 1 I , M„ := | 1 I , M/3 := | B I , Mi := 



For some selfadjoint operator B G L{Hi). The resulting material law is given by 

Mo + 5o-°M« + Oq^M^ + Oq^Mi. 

0^ 
In the situation of Condition 13. 4| we have Po= |000|,Fo= |010| and Qq 
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4 Some Applications 

0\ 

I . Then it is easy to see that Condition 13.41 is satisfied. Furthermore, note that since 
,001/ 

/I 0\ /O 0~ 

one has to compensate the coefficient Mg, one cannot choose Pq = I 010 |,-Fo= I 000 

Vooo/ Vooo 



4 Some Applications 

We want ffi'st to focus on the case Mq = 0. The most simple case of this type is now 

with a = 1 — /3, /3 g]0, 1[, yielding problems of the form 

' d^M^ + Mi+A)u = F, (20) 



for given bounded linear operators M^, Mi, a skew-selfadjoint operator ^ in a Hilbert space 
H, right-hand side F G Hg^Q(R, H) and suitable £» G ]0, oo[. A particular instance of this first 
problem class is the fractional Fokker-Planck equation to be discussed in Section 14.11 A more 
complicated material law including an infinite number of fractional integrals is discussed in 
Section [4.21 within the context of fractional visco-elasticity. Section [4.31 outlines possible ways 
to formulate initial (boundary) value problems for evolutionary equations of the form ()20p . 



4.1 Fractional Reaction-Diffusion Equations 

Let us consider the fractional Fokker-Planck equation, compare e.g. [5]. The operators involved 
can be defined as follows. Let Vt C W^ be an open subset. We denote 

grad,: Coo {^) ^ i^(^) ^ L^{^T 
f ^ (5i/)ig{i,...,„,} 

and div := — grad* and grad := gradj^. The operator matrix 

div' 



A- , 

grad 

is then skew-selfadjoint. Of course other boundary conditions can also be treated, see e.g. 
|11 | [16 } [To]. Let a g]0, 1[. For the material we assume 

V / \fiio iiii 



■^Of course the domain of grad is the Sobolev space Hq{Q,) of weakly differentiable functions satisfying 
generahzed homogeneous Dirichlet boundary conditions. The domain of div is the set of L^(r2)-vector 
fields with distributional divergence in L^(Sl). 
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4.2 An Application to a General Fractional Kelvin-Voigt Model for Visco-Elastic Solids. 

with Ka, yit flu continuous, selfadjoint and strictly positive definite on the respective L^-type 
block component spaces, /ij^ bounded, linear operators, i,k £ {0, 1}. Thus, the block operator 
matrix yielding the (fractional) Fokker-Planck equation looks formally like 

do~"Ka + Hoo div+^oi 
grad + ^10 fJ-u 

By our general theory developed above, the latter operator matrix is continuously invertible 
in the Hilbert space of ffg^Q-functions taking values in L'^{Q,) for sufficiently large ^ > 0. A 
simple row operation yields 

f OQ~°'Ka + fJ'OO - (div +^oi) ^n (gi"ad + ^lo j \ 
\ grad + ^fio /Uii / 

The operator in the first row of the block operator matrix reduces to the following fractional 
convection-diffusion operator 

<9o~"Ka + fiQo - (div +/ioi) fJ'il (g^ad + /iio ' 
= dg'^Ka + {fJ-oo - /Woi^nVio) - div ^j;/ grad - ^oi/^n grad - div/i^/^io- 
The second row represents the flux <^ in terms of the unknown probability density "d: 

^ = -/ij^/gradi? - ^nVioi?- (21) 

Clearly, formally for a = we recover the usual convection-diffusion equation. Speaking in 
terms of heat conduction, which is, up to a choice of units (and interpretation), governed by 
the same type of convection-diffusion equation, (|21|) is a variant of Fourier's law. 



4.2 An Application to a General Fractional Kelvin-Voigt Model for 
Visco-Elastic Solids. 

As another application we consider the material law associated with a fractional Kelvin-Voigt 
model. It is a well-established idea that material laws describing visco-elastic solids may be 
more simple to match with measurements if fractional integrals are admitted, see e.g. [I], [6]. 
The spatial operator 

- Div ' 



^''^ ' Grad 

where Div is the restriction of the tensorial divergence operator div to symmetric tensors of 
order 2 and Grad is the symmetric part of the Jacobian matrix d0u oi u, is skew-selfadjoint if 
suitable boundary conditions at the boundary of the underlying domain C M^ are imposed, 
e.g., similar to the previous example, vanishing of the displacement u on the boundary. This 
choice of boundary condition would amount to replacing Grad by the closure Grad of the 
restriction of Grad to vector fields with smooth components vanishing outside of a compact 
subset of fi. We will not be explicit on the boundary conditions under consideration and 
simply assume that A is considered with a domain such that it becomes skew-selfadjoint. We 
refer to |16l Section 4.2] or |17l Section 2] for more details. The above general perspective to 
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4 Some Applications 

evolutionary equations allow for very general materials to be covered by the approach. 

Let us take a closer look at a fractional Kelvin- Voigt model for visco-elastic material. In this 
model, given an external forcing term /, we have the equation 

doTjdou — DivT = / 

linking the stress tensor field T with the displacement vector field u, where r] is a bounded 
selfadjoint and strictly positive definite operator, accompanied by a material relation of the 
form 

T={C + Dd^)£ (22) 

with £ := Gradti, a € ]0, 1]. Here the case D = would correspond to purely elastic behavior 
and the case a = 1 and D ^ would describe the classical Kelvin- Voigt material. Introducing 
V := Oqu as a new unknown and differentiating (j22p . we obtain 

doijv — DivT = / 
do {C + DdSr^ T = Giadv 

yielding formally an evolutionary equation of the form 

^°^ {C + DdSr^)^[-Giad ))\t)^[o 
Thus, we obtain a system of the form 

where ^ is a skew-selfadjoint realization of the formal block operator matrix given by 

-Div 
- Grad 



and the material law operator 



{C + Dd^y^ 



In order to warrant existence and boundedness of the operator (C -|- DOq) , we impose the 
condition C,D > and 

C + Dq"" > Co (23) 

for some cq G M>o and all sufficiently large g € ]0, oo[. However, in order to show that M (^Oq ) 
satisfies the well-posedness condition, we have to impose a slightly stronger condition. 

Theorem 4.1. Assume that C, D are non-negative selfadjoint operators and such that l^ ^Clq^o, 
lI jjDli^d o-f^ both strictly positive definite, where lq^o and li^d are the canonical embeddings of 
the null space and range of D in the underlying Hilhert space of symmetric 3 x 3-matrices with 
entries in L^ (il), respectively. Then estimate \2'J\) holds and (C -|- D(9q )~ satisfies Condition 
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4.2 An Application to a General Fractional Kelvin-Voigt Model for Visco-Elastic Solids. 



Proof. First of all note that estimate (|23p is a straight forward consequence of Euklid's in- 
equality. Now, we have by decomposition (see also Lemma l3.2p : 

We abbreviate Cjk ■= L*jjCLk,D (j, A; e {0, 1}), Du := lIjjDli^d, 
and W ■.= {^ ~'^oo C'oi \ rpj^^^^ ^^ ^^^ ^j^^^ 



Cio (Cn + i?iia^) ; I Cn + I^iia^ 



With g (9o-^) := (^ao-^^Z^ri^Cn^^n + ij = E:r=o5o-"" (^- V^n ^nV^n ) for Q 
sufficiently large, we have 

Thus, 

{c + Dd^y^ 



^ , 1- 



('-'-"no8„-^e(a„-.)^^^^■•'°■" 







00 /o 

00/0 






/ E ' ('o,-.o)^ (Lry^a,"" (-^cn^vv^ ) '^^ (:?;: 



T I W >° 

11 / V-l.D 
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4 Some Applications 



Hence, with 



V 







^" " ' {io,D ii,D)w(^^^' ^^]W*( '0'^ ^ ' ' 



-1 







^l,D 



and 



/o 



M^ 







V 



( Lo,D Ll,D ) W 



\/i? 



'ii' ( -\/^n Gii\/^n ) V^n ' J ^* l^^D 



for 7 G n := {(1 + n) a;n £ {0, . . . ,\a ^] — 1}} we have 

M (^0-^) = Mo + J^ V^7 + do^M2 (9o"') , 

7en 



where 

oo 

:= ^ ( '-o.D ii,Z) ) W^ 

n=\a~^~\ 







5o"V^n ^o""" 



^iTCiiV^ii V^ 



VF* 



'11 



I* 



Now, for g € M>o sufficiently large we estimate 



\M2{do^)\\ <K, 



do E ^Q-^^o""" -V^n'^iiV^i^i' 



71= I a^ 
oo 



<Ko Y. ^-"+'-""^r = ^0^'-" j^ (^-°i^i)' 



Koe 



1-a 



1-^- 






1 - g-'^Ki 1 - ^~"Ki 






1 - g-°Ki 



1 - ^-"i^i 



< Kog-'' 



K 



[a-n 



1 - g^'^Ki 



^0 (£1 ^ oo) 



where i^r 







Wx D 



'11 



the respective choices 



and K 



Du'CuJD^i 



. Thus, Condition 13.41 is satisfied, for 



Po=(Jo), Fo=([jJ ), andQo = 0. 



D 



Remark 4.2. In |17| a well-posedness result for the fractional Kelvin- Voigt model considered 
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4.3 Initial Boundary Value Problems 
was shown under the condition that a ^ 1/2. 

4.3 Initial Boundary Value Problems 

We now consider initial value problems for equations of the form (|2U|) . Since fractional deriva- 
tives of order /3 G]0, 1[ are non-local, i.e., we have memory effects, prescribing a pre-history 
appears to be more appropriate. This amounts to imposing Hq^q (M, //)-data prior to the 
initial time 0, i.e., "the pre-history", and is clearly covered by the above. 

However, if initial conditions are deemed appropriate, then there appear to be two possible 
ways of implementing them. The first option would be to consider with F G Hg^Q{M,H), 
W^ e M„ [H] , 

d^MaU + MiU + AU = F + 6^Wa (24) 

as a straight-forward generalization of the case /? = 1, see |101 Theorem 6.2.9]. This translates 
to 

dodQ^'MaU + MiU + AU = F + S^Wa, 

which following the rationale in |101 Proof of Theorem 6.2.9] would be in point- wise terms an 
implementation of the initial jump condition 

(ao-"M„c/ - x«^„ ® w^) (0+) - (ao-"M„f/ - x^^, ® w^) (o-) = o, 

where the limits are taken in H^i := H^i{A + 1), the completion of H with respect to the 
norm |(^ + 1)~^ • \h, see also Section [2.21 If F = on IR<o then causality yields 



^do'^M^U - x^^, ^ Wa) (0-) = in H_i 
and so, in this case an initial condition of the form 

[d^^M^U) (0+) = Wo, in H^i 

results. 

Alternatively, let 

Va € M, [H] 

and G € i^g,o (^^ H) be given and consider 



5^ [M^V - x^^, ^ Va) + MiV + AV = G. (25) 

For the case a = 0, f3 = l,G = F and Vq = Wq this would be a re- formulation of (|24|) . This 
justifies taking ()25p as an appropriate alternative generalization. Problem ([25]) is equivalent 
to the evolutionary problem 

d^M^V + MiV + AV = G + 5o^Xh>o ^ ^"' 

= G + ao"°(5 0Kv. 

Both approaches are clearly closely related. Indeed, for G = d^^F and Va = Wa we see by 
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